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ABSTRACT 

We present the results of a multisite photometric campaign on the pulsating sdB star Bal- 
loon 090100001. The star is one of the two known hybrid hot subdwarfs with both long- and 
short-period oscillations, theoretically attributed to g- and p-modes. The campaign involved 
eight telescopes with three obtaining UBVR data, four B-band data, and one Stromgren uvhy 
photometry. The campaign covered 48 nights, providing a temporal resolution of 0.36 /iHz 
with a detection threshold of about 0.2 mmag in S-filter data. 

Balloon 090100001 has the richest pulsation spectrum of any known pulsating subdwarf 
B star and our analysis detected 1 14 frequencies including 97 independent and 17 combination 
ones. Most of the 24 g-mode frequencies are between 0.1 and 0.4 mHz. Of the remaining 73, 
presumably p-modes, 72 group into four distinct regions near 2.8, 3.8, 4.7 and 5.5 mHz. The 
density of frequencies requires that some modes must have degrees £ larger than 2. The modes 
in the 2.8 mHz region have the largest amplitudes. The strongest mode (/i) is most likely 
radial while the remaining ones in this region form two nearly symmetric multiplets: a triplet 
and quintuplet, attributed to rotationally split t = \ and 2 modes, respectively. We find clear 
increases of splitting in both multiplets between the 2004 and 2005 observing campaigns, 
amounting to ^15% on average. The observed splittings imply that the rotational rate in Bal09 
depends on stellar latitude and is the fastest on the equator. We also speculate on the possible 
reasons for the changes of splitting. The only plausible explanation we find are torsional 
oscillations. This hypothesis, however, need to be verified in the future by detailed modelling. 
In this context, it is very important to monitor the splittings on a longer time scale as their 
behaviour may help to explain this interesting phenomenon. 

The amplitudes of almost all terms detected both in 2004 and 2005 were found to vary. 
This is evident even during one season; for example, amplitudes of modes /§ and /c were 
found to change by a factor of 2-3 within about 50 days during 2005. 

We use a small grid of models to constrain the main mode iji), which most likely rep- 
resents the radial fundamental pulsation. The groups of p-mode frequencies appear to lie in 
the vicinity of consecutive radial overtones, up to the third one. Despite the large number of 
g-mode frequencies observed, we failed to identify them, most likely because of the disrup- 
tion of asymptotic behaviour by mode trapping. The observed frequencies were not, however, 
fully exploited in terms of seismic analysis which should be done in the future with a larger 
grid of reUable evolutionary models of hot subdwarfs. 

Key words: oscillations - subdwarf - stars: individual: Balloon 090100001. 
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1 INTRODUCTION 

lKilkennvet"al] ( ll997l) discovered the first group of B-type hot sub- 
dwarf (sdB) pulsators, which they termed EC 14026 after the pro- 
totype and are now designated as V361 Hya stars, and referred to 
as sdBV stars. This discovery opened the possibility of constrain- 
ing the interiors of sdB stars through asteroseismology, which, in 
turn, may help to better understand their still controversial evolu- 
tion. So far, about 40 class members have been found, with periods 
ranging from 60 to 600 s and typical amplitudes of 10 mmag, with 
some exce ptions of higher amplit udes up to 60 mmag. Theoretical 
studies by jcharpinet et al] ( Il99'/) demonstrated that their oscilla- 
tions could be explained by acoustic (p-) modes of low degree, I, 
and low radial orders, n, driven by an opacity bump caused by a 
local enhancement of iron-group elements. 

In 2002 ano ther class of pulsating sdB stars was found 
dOreen et al. I l2003h . The new class is characterised by lower am- 
plitudes (1 -2 mmag) and longer periods (~1 hour). The new vari- 
able stars are referre d to as long-p e riod s dBs, or PG1716 stars, 
after that prototype. iFontaine et al.l j2003l) determined that their 
variability could be explained in terms of high radial order grav- 
ity {g-) modes. The same n mechanism that operates in V 361 Hya 
stars also works in PG1716 stars, although theoretical unstable 
modes were at slightl y lower effective temperatures than observed. 
llefferv&Saiol ( l2006l) , using OP opacity tables and considering an 



Balloon 0901 00001 : 2005 campaign 3 



enhancement of both iron and nickel in the driving zone, were able 
to obtain unstable modes at the same effective temperatures as the 
observed ones. Additionally, their non-adiabatic analysis led to si- 
multaneous instability of p- and (;-modes, which can account for 
the two known hybrid sdB pulsators. 

The two classes (V361 Hya and PG1716 stars) of pulsat- 
ing sdBs have slightly different effective temperatures and surface 
gravities, the former being hotter and denser than the latter. Their 
pulsating behaviour resembles two other pairs of pulsating stars 
showing p- and g-modes: /3 Cephei and Slowly Pulsating B stars 
and 5 Scuti and 7 Doradus stars. 



2 BALLOON 090100001 

Balloon 090100001 (hereafter Bal09) is one of the brightest pul- 
sating sdB stars discovere d to date. It was found to be an sdB 
star by [ Bixler et ^ jl99lh . and a short-period pulsating sdB by 
lOreiro et al 1 2004 ). It has relatively long periods and large am- 
plitudes of oscillation, which makes this object an excellent can- 
didate for follow-up photometry. Bal09 was observed in a long- 
time base campaign independently by two of us, in August and 
September 2004 using the 60-cm telescope at Mt. Suhora Observa- 
tory (AB) and the 80-cm IAC80 telescope at Tenerife (RO). Four 
wide-band filters, UBVR, were used in the former jBaran et"ail 
I2005L her eafter BarOS), whil e only B-filter data were obtained in 
the latter ( lOreiro et alj|2005l , hereafter Ore05). The use of a com- 
mon filter allowed an analysis of the combined data, from which al- 
most 50 frequencies were detected (Baran et al. 2006). These cam- 
paigns revealed a rich spectrum of frequencies clustering in groups 
within the p-mode domain, and having lower amplitudes as the fre- 
quency increases. A high-amplitude frequency (~ 60 mmag) dom- 
inates the spectrum, with an equally-spaced triplet nearby. Assum- 
ing that the triplet splitting is caused by rotation, a rotational pe- 
riod near 7 days was derived. Interestingly, variability near 4 mmag 
was also detected in the low-frequency region, which is typical for 
PG1716 stars. Thus, Bal09 turned out to be a hybrid object and 
became very interesting from the seismological point of view. The 
star was the second hybrid sdB pulsa tor discovered, the first be- 
ing HS 0702-1-6043 l iSchuh et alj|2006l) , which shows a single fre- 
quency in the g-mode region. 
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3 2005 PHOTOMETRIC CAMPAIGN 

Despite the frequency resolution achieved with the 2004 observa- 
tions, the spectrum of Bal 09 suffered from strong aliasing, hinder- 
ing a clear frequency identification at lower amplitudes. Closely- 
spaced frequencies and/or amplitude variability made the situation 
even worse. 

In order to get through these problems, a multisite campaign 
was organized in the summer of 2005. It involved eight tele- 
scopes at different longitudes, and spanned almost two months 
(10th August -27th September). All sites but one used CCD cam- 
eras with Johnson B filter (or very similar), although some of 
them performed multicolour photometry in UBVR filters. More- 
over, Stromgren uvby photometry was acquired with a photoelec- 
tric photometer attached to the 90-cm telescope in Sierra Nevada 
Observatory. The log of observations, including sites, observers and 
other relevant information is given in Table [T] 
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Figure 1. Observational data of Bal 09 in B and the average of b and v 
filters obtained in August and September 2005. Each panel covers one day 
and the numbers in the bottom left or right comers represent the three last 
integer digits of HJD for that panel. Note that consecutive days go from top 
to bottom. The range in ordinate amounts to 240 mma. 



3.1 Data reduction 

Raw CCD images were calibrated in a standard way by removing 
instrumental effects using bias, dark and flat field images. Mag- 
nitudes of stars were extracted using aperture photometry, but by 
means of different reduction pipelines: data collected at I AC 80 
(Tene rife) were reduc ed with the Real Time Photometry (RTP) soft- 
ware (0stensen 2000); IDL photometry packages were applied to 
the Lulin (Taiwan) data; photometry from Sobaeksan (Korea) and 
Mt. Lemmon (Arizona) were proccessed with IRAF routines, while 
other CCD data were reduced by means of the DAOPHOT package 
with a DAGGROW supplement l lStetson| [T987, 1990). Data from 
the four-channel photoelectric photometer attached to the 90-cm 
telescope in Sierra Nevada Observatory were corrected with nightly 
extinction coefficients derived from a nearby star. 

Differential photometry for all filters was obtained using the 
reference star GSC 02248.00063, the most suitable due to its lo- 
cation in the field of view and magnitude. However, as it is much 
cooler than Bal 09, all data were corrected for second-order extinc- 
tion effects. If, after this step, the light curve still showed long- 
period variations (on time scales much longer than a few hours), a 
spline fit was applied and subtracted from the data to remove those 
trends. The resulting light curves, obtained for each site separately, 
were then shifted to the same mean, converted to relative flux and 
combined. A similar conversion to relative fluxes has been done for 
the 2004 data. Consequently, all amplitudes provided in this paper 
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Table 1. Details on the 2005 campaign. UH stands for University of Hawaii at Hilo. Other observer's abbreviations correspond to initials of the co-authors. 



Telescope 


r liter 


rl.l I J Slail 


rl.l IJ enu 


XT- 
Nights 


Hours 


Observers 


lAC Teide, 0.8 m 


U 


620.37 


628.75 


7 


44.2 


RO, JMGP. TM, LF 




B 


593.54 


628.75 


18 


115.2 






V 


620.37 


628.75 


7 


44.2 






R 


620.37 


628.75 


7 


44.1 




Mt. Suhora, 0.6 m 


U 


594.47 


639.51 


34 


226.9 


AB, MS, DK, SZ 




]j 


594 47 


639 51 


36 


216.6 






V 


594.47 


639.51 


34 


226.9 






R 


594.47 


639.51 


34 


226.9 




Baker, 0.4 m 


U 


594.68 


624.91 


11 


79.2 


MR, AYZ, SH, JRE 




B 


594.68 


625.87 


13 


92.4 






V 


594.68 


594.79 


1 


2.7 






R 


594.68 


624.96 


11 


80.0 




Lulin, 1 m 


B 


609.01 


613.13 


5 


21.9 


WPC, HTL, FYH 


Sobaeksan, 0.6 m 


B 


618.01 


630.03 


5 


23.3 


SLK, KJC 


Mauna Kea, 0.6 m 


B 


621.87 


631.02 


7 


36.7 


AB, RC, UH students 


Mt. Lemmon, 1 m 


B 


636.84 


640.63 


5 


29.0 


SLK 


Sierra Nevada, 0.9 m 


uvby 


605.41 


618.62 


14 


105.6 


CRL, RG 



are given in the units of 1/1000 of the relative flux (mma). In Fig.[T] 
the entire B and bv averaged hght curves are presented. 



3.2 Re-analysis of the combined 2004 data 

As previously described, two independent photometric data sets of 
Bal09 were obtained in 2004. As our analysis of the 2005 cam- 
paign data revealed amplitude variations (see Section [33] l. we de- 
cided to combine the 2004 B-filter data and re-analyse them in the 
same way as the 2005 campaign data. The 2004 data of Bar05 are 
better distributed in time and cover a longer time interval than those 
of Ore05, allowing for a better distinction of close frequencies. On 
the other hand, the data of Ore05 have better precision and were 
carried out with shorter time sampling. In order to balance the ac- 
curacy and resolution in the combined data, we decided to average 
three consecutive datapoints in the data of Ore05 before combining 
them together. 

Section [331 explains details of the analysis of the 2005 data. 
For the sake of consistency, the combined 2004 data were analysed 
in the same way. The new feature was the inclusion of the rates 
of amplitude changes in the non-linear least squares used to sub- 
tract periodic terms found by means of Fourier periodogram. The 
rates of these changes as well as the amplitudes and phases of all 
frequencies detected in the 2004 data are reported in Table|2] In to- 
tal, we detected 56 frequencies of which 41 (12 in g-mode and 29 
inp-mode regions) are independent and 15 are combination terms. 
In comparison to the analysis of Bar05, we found 26 more frequen- 
cies. This is mainly the consequence of lowering the detection level 
in the combined data. Of the 30 modes given by Bar05, we did not 
detect /e and /h in the g-mode region. This creates uncertainty 
as to whether they are intrinsic to Bal09, but it is also possible 
that their amplitudes were changing quite fast, making them unde - 
tectable in the combined data (see for example iReed et al]|2007ah . 
This is only possible if their amplitudes were very low at the begin- 
ning of the 2004 run when most of the Ore05 data were acquired. 
It is worth noting that the combined data revealed a fourth region 



of p-mode frequencies around 5.5 mHz. Frequencies in this region 
were reported by Ore05, but not detected by Bar05. 

In Table|2]and throughout this paper, the notation of frequen- 
cies introduced by Bar05 is used and extended for new modes: fre- 
quencies in the g-mode region are annotated by letters while num- 
bers are assigned to frequencies in the p-mode region. When a fre- 
quency was recovered at an alias frequency it is marked with '+' or 
' — '; for example /12+ = fi2 + 1 d"^ and /12- = fi2 — 1 d"^. Ta- 
ble 2 of Bar05 indicated two possible aliases for a single frequency 
designated as /ig. We shall denote the two possibilities as /ig/ and 
/18" , the frequency of the former being lower than the latter. Ac- 
cording to this notation, the frequency listed as /14 in Table 2 of 
iBaran et al.. C2008.) should be designated as /ig// . 



3.3 Fourier analysis of the 2005 campaign data 

Both the 2004 and 2005 data for Bal09 cover similar time intervals: 
nearly 40 days in 2004 and 50 days in 2005. The resulting resolu- 
tion is therefore similar, but the improved duty cycle of the 2005 
data results in a lower detection threshold and improved aliasing. 
In terms of time coverage, the 2005 UBVR and Stromgren uvby 
data are complementary. Therefore, we decided to combine them. 
In order to minimize the wavelength dependence of pulsation am- 
plitudes, we combined the U data with the mean of the u and v 
(hereafter Uuv) data, the B data with the average of the v and b 
(hereafter Bvb) data, and the V and y (hereafter Vy) data. The 
Bvb data are the most numerous of the four data sets and so pro- 
vide the best spectral window which is shown in Fig.|2] 

In Fig. [3] the Fourier spectrum for the Bvb data is shown. 
Since oscillations in Bal09 are dominated by frequencies at 
2.8 mHz, especially the strongest one with an amplitude near 
50 mma, for clarity Fig. |3] presents the amplitude spectrum af- 
ter removing the three highest-amplitude frequencies (/i, /2 and 
/a). It can be seen that the same groups of frequencies are present 
that were in the 2004 campaign (Bar05, Ore05): the low-frequency 
region (mainly around 0.3 mHz), that of the dominant frequency 
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Table 2. A list of modes detected in the combined 2004 B filter data, while allowing for linear amplitude changes. Amplitudes and phases are given for HJD 
= 2453 246. Numbers in brackets indicate the formal rm.s. errors. 



# Frequency Ampl. AAlAt Phase # Frequency Ampl. dA/dt Phase 

[mHz] A [mma] [mma/d] [rad] [mHz] A [mma] [mma/d] [rad] 



fl 


0.229566(24) 


0.57 


+0.013 


0.12(17) 


/26 


3.822911(21) 


0.70 


-0.002 


0.14(13) 


/d 


0.239973(07) 


1.96 


+0.010 


5.80(05) 












/f 


0.246315(22) 


0.61 


+0.003 


2.64(15) 


fu- 


4.633620(25) 


0.57 


0.000 


1.26(17) 


/a 


0.272379(05) 


2.75 


-0.028 


5.66(03) 


/21 


4.644467(19) 


1.14 


-0.034 


1.96(10) 


/g 


0.298886(19) 


0.67 


+0.003 


6.23(13) 


/27 


4.651166(25) 


0.65 


-0.005 


2.20(15) 


fc 


0.325671(12) 


1.25 


+0.052 


1.37(09) 


/28 


4.655560(28) 


0.80 


-0.022 


4.35(15) 


h 


0.331210(16) 


0.81 


+0.007 


1.42(11) 


/29 


4.660107(22) 


0.76 


-0.009 


2.70(13) 


/b 


0.365807(05) 


2.66 


-0.023 


1.98(03) 


/l3 


4.661418(14) 


1.07 


+0.005 


2.04(09) 


/k 


0.397189(26) 


0.49 


+0.011 


0.48(20) 


/30 


4.670320(17) 


0.92 


-0.002 


4.24(11) 


/l 


0.631078(29) 


0.44 


+0.011 


0.94(22) 


/31 


4.676034(34) 


0.65 


—0.026 


2.57(21) 


/m 


0.684352(21) 


0.60 


+0.002 


4.77(14) 


/12+ 


4.680947(21) 


0.78 


-0.006 


4.00(12) 


/n 


0.833085(25) 


0.53 


-0.006 


4.84(17) 






















/32 


5.506010(28) 


0.54 


-0.013 


1.16(18) 


h-h 


0.968616 


0.43 


+0.001 


5.21(20) 


/33 


5.532958(21) 


0.61 


+0.009 


3.86(15) 












/34 


5.553408(27) 


0.48 


-0.002 


0.11(18) 


fl - /b 


2.441662 


0.80 


+0.001 


1.11(11) 


/35 


5.605820(28) 


0.47 


+0.002 


4.13(19) 


fl - fc 


2.481799 


0.43 


+0.016 


2.13(19) 












fl - /a 


2.535091 


0.41 


+0.009 


3.57(20) 


2/1 


5.6149394 


5.67 


+0.003 


3.874(16) 












/1+/2 


5.6307087 


4.33 


0.000 


1.793(21) 


fl 


2.8074697(02) 


52.83 


+0.058 


2.6938(17) 


/1+/3 


5.6322777 


2.61 


-0.022 


0.94(04) 


f2 


2.8232390(06) 


20.39 


+0.032 


0.602(05) 


/1+/4 


5.6338393 


0.99 


-0.003 


4.10(10) 


/3 


2.8248080(11) 


12.13 


-0.055 


5.997(07) 


2/2 


5.6464780 


0.80 


+0.001 


5.75(12) 


/4 


2.8263696(27) 


5.09 


-0.034 


2.784(18) 


/2 + /3 


5.6480470 


0.66 


-0.002 


4.87(14) 


/5 


2.853958(08) 


1.73 


-0.001 


5.19(5) 


/1+/5 


5.661428 


0.53 


+0.002 


0.44(17) 


/7 


2.855710(10) 


1.38 


+0.009 


2.27(7) 












/6 


2.858534(14) 


1.15 


-0.002 


4.86(9) 


/l +/8 


6.583556 


0.55 


-0.013 


1.90(15) 


fl5 


3.763702(25) 


0.61 


+0.009 


1.51(17) 


3/1 


8.4224091 


0.71 


-0.010 


5.37(12) 


f23 


3.764802(30) 


0.64 


-0.003 


2.54(16) 


2/1 + /2 


8.4381784 


0.49 


+0.005 


2.83(18) 


fs 


3.776086(05) 


4.50 


-0.045 


0.914(24) 


2/1 + fs 


8.4397475 


0.40 


-0.009 


2.05(21) 


/9 


3.786724(13) 


1.33 


-0.007 


2.51(08) 












/24 


3.791843(18) 


1.00 


-0.028 


3.61(11) 






Nobs 


12115 




/lO 


3.795572(13) 


1.05 


+0.003 


1.60(09) 




Detection threshold [ 


mma] 


0.42 




/l8' 


3.797813(23) 


0.68 


+0.007 


4.94(16) 




(T4, cr^A/dt [mma, mma/d] 


0.09 


0.007 


/25 


3.821887(28) 


0.54 


+0.014 


1.22(20) 




Residual S.Dev. [ 


mma] 


6.46 






2.8 2.84 
Frequency [mHz] 

Figure 2. Spectral window of 2005 Bvb data of Bal 09 shifted to the fre- 
quency of the strongest mode. 



(at 2.8 mHz) and three otlier groups (3.8, 4.7, 5.5 mHz) separated 
roughly by the same frequency distance (0.8-1.0 mHz). However, 



their amplitudes have changed; the low-frequency pulsations, in 
particular, increased their amplitudes between 2004 and 2005. 

^From the 2004 data Bar05 found only a weak evidence for 
amplitude variations, but those data were only sensitive to varia- 
tions on a time scale of one or two weeks. However, amplitude vari- 
ations are obvious between the 2004 and 2005 data sets. In order to 
infer the character of amplitude changes, we divided the 2005 Bvb 
data into 33 subsets. The subsets contained data covering time inter- 
vals of 10 days in order to ensure that the components of multiplets 
at 2.8 mHz will be resolved. The subsets, however, largely overlap. 
The first subset started at the beginning of Bvb data, the next one 
a night later than the first one and so on. The subsets were anal- 
ysed independently and the amplitudes for only the six strongest 
frequencies (/i through /4, /g, and /c) were used. They are plot- 
ted in Fig. |4l Of the six frequecies shown in Fig. |4l only shows 
no clear evidence for amplitude variations. 

All four {Uvb, Bvb, Vy and R) data sets were analysed in- 
dependently. The procedure of extraction of periodic terms was a 
typical pre-whitening procedure, consisting of the following steps: 
(i) finding the frequency of the most significant peak in the Fourier 
periodogram calculated using the residuals from previous fit (or 
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Frequency [mHz] 



Figure 3. Fourier spectram of the 2004 B (top) and 2005 Bvb data (bottom) 
after removing the three strongest frequencies near 2.8 mHz (/i, /2 and 
/3, see Tables|2]and[3). Regions attributed to p-modes, g-modes and their 
combinations are labeled with 'p', 'g', and 'c', respectively. 



the original light curve when starting the procedure), (ii) fitting a 
truncated Fourier series by means of least-squares with frequencies 
fixed; all previously found terms, also the combination ones, were 
included in the fit, (iii) improving the frequencies, phases and am- 
plitudes of these terms using non-linear least-squares fit. The resid- 
uals from such a fit were used in step (i). The steps (i)-(iii) were 
repeated until no significant peaks in the Fourier periodogram were 
found. 

In view of significant amplitude variations observed in Bal 09, 
in the final analysis we decided to allow for amplitude changes, at 
least for the strongest modes. Therefore, for the strongest terms the 
rates of amplitude changes were added in the step (iii) as free pa- 
rameters, with A-AIAt = as starting values. By the strongest terms 
we mean those which had amplitudes at least ~3.5 times larger 
than the detection threshold in a given data set, i.e. 2.0, 0.7, 1.6 
and 1.5 mma, for Uuv, Bvb, Vy and R data, respectively. (If for 
a given term the rate of amplitude change was derived, the result- 
ing value is reported in Table [3]) For the remaining terms, dA/dt 
was fixed to zero. Fixing the rates to zero for the small-amplitude 
terms allowed us to avoid problems with convergence which oth- 
erwise occurred. For Uuv, Vy and R, the process was continued 
until all terms with signal-to-noise (S/N) above 4 were extracted. 
The resulting frequencies, phases, amplitudes and rates of ampli- 
tude change are given in Table [3] 



E 
E_ 

0) 
T3 
13 

■q. 
E 
< 



— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 



-H 1 1 1 H 



-H 1 1 1 H 



595 600 605 610 615 620 625 630 635 640 
HJD - 2453000.0 

Figure 4. Amplitudes of the six strongest frequencies in Bal 09 for the 2005 
Bvb data. ^From top to bottom the frequencies are /i, /2, /a, fi, fs and 
fc (see Table|3). For f-j, and /c the dashed lines con'espond to changes 
calculated from the derived values of A and dA/dt (Table |3), for the re- 
maining three frequencies they come from a more complicated model (see 
text for an explanation). 



For the Bvb data, we used the same detection limit (S/N > 
4), but the process required a modification as the non-linear least- 
squares routines did not reach convergence once a large number of 
frequencies were fit. Therefore, we first obtained a preliminary so- 
lution which included all frequencies with amplitudes higher than 
1 mma. Next, we took advantage of the fact that in Bal 09 frequen- 
cies appear in groups separated by several hundred fiHz. Thus, 
we could safely assume that during least-squares fits, frequencies 
from one group do not affect the solution for frequencies of the 
other groups. We obtained the final solution for Bvb data using 
six groups. In this process, we first subtracted frequencies from 



the other groups taken from the preliminary solution and used the 
residuals as a starting point for the extraction of frequencies in the 
considered group. 

As can be seen from Fig. [4] the amplitude variations can be 
more complicated than linear. We may therefore expect the oc- 
curence of residual signals in the vicinity of strong terms result- 
ing from amplitude changes unaccounted for by the model we use, 
i.e., linear amplitude change. We detected several such terms (we 
call them 'ghost' terms) near /i, /2 and /g. Since they are arte- 
facts caused by amplitude changes not accounted for by a constant 
dA/dt, they are not listed in Table |3] Having their amplitudes, fre- 
quencies and phases, we can, however, combine them with the main 
term to show the apparent amplitude changes. This is shown in 
Fig. E] for /i, /2 and fs- For the remaining terms shown in this 
figure the 'ghost' terms were not detected which means that the 
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Table 3. Modes detected in 2005 data in the Uuv, Bvb, Vy and R data. Numbers in parentheses indicate the r.m.s. en'ors. Phases and amplitudes are given 
for epoch HJD 2 453 615. 



# 


Mean freq. 




Amplitude, A [mma 


,dA/dt [mma/d] 






Phase [rad] 






[mHz] 


IJuv 


Bvb 


Vy 


R 


Uuv 


Bvb 


Vy 


R 


/o 


0.096902(09) 




0.70,-0.031 








4.86(07) 






/p 


0.139161(22) 






0.77, — 








1.15(13) 




/q 


0.142459(10) 


... 


0.59, — 








3.37(07) 






/r 


0.144996(22) 


... 




0.79, — 








5.65(13) 




h 


0.169849(10) 




0.58,— 








4.09(07) 






h 


0.175656(08) 




0.68,— 


0.88, — 






0.90(06) 


0.56(12) 




fu 


0.192869(21) 


... 




0.83, — 








3.41(12) 




fv 


0.198935(11) 


... 


0.54, — 








5.89(08) 






h 


0.229549(09) 




0.85,-0.027 




0.49, — 




5.01(06) 




4.95(19) 


/d 


0.239971(07) 


1.68, — 


1.08,+0.019 


0.99, — 


0.85,— 


5.81(08) 


5.99(04) 


5.98(10) 


5.73(11) 


J r 


0.246301(05) 


1.80, 


1.52,+0.018 


1.66,-0.004 


1.37^ 


1.79(07) 


1.613(27) 


1.58(06) 


1.57(07) 


/a 


0.272463(05) 


1.67,— 


1.72,-0.002 


1.45,— 


1.49,— 


1.37(08) 


1.473(23) 


1.58(07) 


1.57(07) 


fw 


0.29380(03) 








0.52, — 








2.20(18) 


/g+ 


0.310575(11) 


0.85,— 


0.72,-0.008 


0.56, — 


0.58,— 


4.57(15) 


4.66(06) 


4.70(18) 


4.88(16) 


fc 


0.3256089(18) 


4.84,-0.077 


4.44,-0.069 


3.82,-0.019 


3.71,-0.053 


4.590(27) 


4.609(10) 


4.609(27) 


4.613(26) 


h 


0.331182(10) 


0.71, — 


0.75,-0.002 


0.79, — 


0.61, — 


0.87(18) 


6.19(05) 


0.58(13) 


6.25(15) 


/b 


0.365805(03) 


3.36,-0.014 


2.65,-0.013 


2.52,-0.006 


2.64,+0.002 


6.20(04) 


6.180(15) 


0.06(04) 


6.22(04) 


/k 


0.397232(09) 


0.67, — 


0.84, 0.000 


0.97, — 


0.56, — 


3.01(19) 




2.63(1 1) 


2.52(17) 


f, 

JL 


63074nfl 51 




0.42 


0.49^ 






1.36(09) 


1.64(21) 




/m 


0.684405(14) 




0.47, — 


0.56, — 


0.38,— 




5.02(08) 


4.90(18) 


4.91(25) 


/n 


0.833090(14) 


0.59,— 


0.43, — 






4.67(22) 


4.77(09) 






/x 


1.845938(27) 




0.20, — 








4.48(19) 






/y 


1.95968(04) 


0.59,— 








0.55(22) 








/z 


2.296144(22) 




0.25,— 








4.23(15) 






/i - /k 


2.410232 




0.22, — 








3.68(18) 






/i -/b 


2.441659 


0.86, — 


0.69, — 


0.51,— 


0.47, — 


6.06(14) 


0.14(06) 


6.23(18) 


6.14(18) 


/i - /c 


2.4818552 


0.81,— 


0.77,-0.013 


0.93,— 


0.63,— 


1.81(14) 




1 74(10) 


1.78(14) 


f ) — fr. 
/3 /C 






20 








4.75(20) 






J AA 


2 7n7505f771 




0.20 








3.34(19) 






2/1 - /3 


2.7900921 




0.20. — 








1.31(19) 






fl 


2.80746412(16) 


59.89,-0.048 


45.48,-0.014 


41.26,+0.019 


40.08,-0.019 


5.912(02) 


5.934(01) 


5.927(03) 


5.934(03) 


h 


2.8230026(08) 


11.16,-0.013 


9.49,-0.033 


8.61,-0.012 


8.27,-0.019 


0.638(12) 


0.658(05) 


0.653(12) 


0.626(12) 


h 


2.8248361(06) 


17.33,+0.029 


13.93,+0.041 


13.49,+0.015 


12.88,+0.042 


2.441(08) 


2.457(03) 


2.454(08) 


2.428(08) 


U 


2.8265995(22) 


4.11, +0.035 


3.54,+0.003 


3.16,+0.001 


2.99,+0.005 


2.16(03) 


2.095(12) 


2.01(04) 


2.11(03) 


h 


2.853383(09) 


1.06,— 


0.92,-0.011 


0.58,— 


0.54, — 


1.00(12) 


0.95(05) 


0.61(18) 


1.13(18) 


fr 


2.855571(10) 


0.66, — 


0.85,-0.008 


0.82, — 


0.68,— 


5.16(20) 


4.99(05) 


5.00(13) 


4.91(14) 




2.857126(13) 


1.05,— 


0.59, — 


0.69, — 


0.59, — 


0.33(13) 


0.54(07) 


0.73(15) 


0.80(17) 


k 


2.858725(10) 


0.91,— 


0.88,-0.001 


0.76, — 


0.96. — 


5.87(14) 


5.74(05) 


5.87(14) 


5.67(10) 


hi 


2.860737(16) 




0.38,— 








4.94(11) 






fss 


3.036271(21) 




0.28,— 








4.96(13) 






fi + fc 


3.1330730 




0.22, — 








0.51(18) 






/i +/b 


3.173269 




0.18.— 








2.36(22) 







model with constant rate of amplitude change is good enough. In- 
deed, there are no large discrepancies from linear amplitude change 
for /s, /4 and fc- 

The results of frequency extraction are presented in the fol- 
lowing subsections and the full list is given in Table [S] The largest 
number of frequencies, 103, was detected in the Bvb data, as ex- 
pected. In each of the remaining three datasets over 30 frequencies 
were detected. The total number of frequencies presented in Table 
[3]is 1 14, since some of them which were not detected in Bvb data 
were found in other datasets. Of these 1 14 frequencies, 17 are com- 
bination terms while the remaining 97 may represent independent 
modes. 



As indicated by Fig. [4] the amplitudes of some terms change 
by a factor of two or more during 50 days (e.g. fs and .fc). As 
such, we cannot exclude the possibility that some nearby frequen- 
cies given in Table [3] are not independent but are artefacts caused 
by amplitude variations. Another problem is aliasing. This is par- 
ticularly bad near 3.8, 4.7 and 5.5 mHz where low-amplitude fre- 
quencies are near to high-amplitude ones and the frequency density 
is very high. Again, it is possible that some of these frequencies 
are in error by 1 (sidereal day)~^. Furthermore, some frequencies 
may also have changed. While these changes probably did not af- 
fect the analysis of any given year data, season-to-season changes 
are evident for some components of rotationally split multiplets. 
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Table 3. continued. 



# 


Mean freq. 


Amplitude, A [mma], AAI dt [mma/d] 






Phase [rad] 






[mHz] 


Uuv Bvb Vy 


R 


U uv 


Bvb Vy 


R 



/39 


7741 QOM 


g4 


49 


64 4 8 1 n ^'l 




3 92(16) 




/8 




9 '^0 —0 006 


9 no ^0 09'^ 


9 96 —0 016 1 Q8 -1-0 008 4. OdiOt^^ 


4 03(03) 


^ Q'^f'O'SI 
J .yj\\j~j ) 


4 09(05'! 


/40 


T 77748nM(S^ 




0.46^ 




1.49(10) 






J41 








0.52 














0.37 




2 71(111 














0.55 






2.25(18) 


/44 = /g? 


3.786291(17) 




0.48, — 


0.43,— 


4.14(09) 




4.21(23) 


/45 = /g? 


3.787153(21) 




0.37, — 




4.18(14) 






/46 


3.791291(15) 




0.58, — 


0.75,— 


2.70(08) 




2.47(14) 


/24 


3.791581(12) 


1.28,— 


0.78,+0.030 


0.74,— 0.57,— 0.91(10) 


0.76(08) 


1.10(14) 


0.56(19) 


/47 


3.792387(22) 




0.36, — 




0.62(13) 






/48 


3.793419(20) 




0.35,— 




4.28(12) 






/49 


3.796825(20) 


0.61,— 


0.39, — 


... 5.45(21) 


5.79(12) 






/19 


3.80520(04) 






0.51,— 




1.09(20) 




/20 


3.806574(22) 




0.30, — 




0.86(14) 






/50 


3.807605(19) 




0.36, — 




5.11(12) 






/18" 


3.808945(15) 




0.44, — 




4.45(09) 






/5I 


3.812194(17) 




0.37, — 




0.97(11) 






/52 = /26? 


3.822485(17) 




0.34, — 




5.13(12) 






/53 


3.824759(12) 




0.55, — 


0.48,— 


1.71(07) 




1.92(20) 


hi 


3.829725(16) 




0.42, — 




2.17(12) 






/55 


3.840824(23) 




0.28, — 




3.15(17) 






fee 
J 56 


4 636265f26i 




0.25 




3.45(17) 






J57 






42 




3 87(10) 






J58 






36 


42 


3 37(12) 




3 87(22) 


fen 






47 


59 


3 17(10) 


2 93(17) 




fan 


d f^d7dl^(^ f,\ 




49 




2 91(10) 






/61 — 711 ■ 


M-.O'-fJl iJ\Z.o) 




27 




1 SQn7'i 






/62 — /l 1 






33 










/63 


4 646299f1 6i 




0.54^ 




4.73(08) 






J 64 


4.647062(28) 




0.28, 




1.77(18) 








4.649385(22) 


0.64, — 


0.29, — 


... 0.11(20) 


0.16(15) 






fee 


4.65294(04) 


0.66, — 




... 0.72(20) 








/67 


4.653491(18) 




0.44, — 




4.11(12) 






/68 


4.654759(20) 


0.74, — 


0.34, — 


... 4.14(18) 


3.71(13) 






/12- 


4.658007(19) 




0.41,— 




4.32(12) 






/69 = /l3? 


4.66097(03) 






0.53, — 




4.14(19) 




/to = /l3? 


4.661708(15) 




0.40, — 




2.02(10) 






/71 


4.663070(25) 






0.62,— 






3.36(15) 


/72 


4.665422(15) 




0.55,— 


0.46, — 


4.36(08) 


4.58(23) 




/73 


4.666868(15) 


0.64, — 


0.58,— 


0.37, — 0.55, — 5.96(20) 


6.01(07) 


6.10(28) 


6.25(17) 


/74 


4.668308(18) 




0.37, — 




4.22(11) 






/75 


4.671570(23) 




0.29, — 




0.75(15) 






/76 


4.674040(26) 




0.24, — 




1.98(17) 






/77 


4.683263(23) 




0.28, — 




0.57(15) 






/78 


4.708165(26) 




0.21,— 




0.25(18) 







as discussed in Section [53] The changes of frequencies between 
2004 and 2005 could be as large as 0.6 /iHz (see Table |4j, which 
may pose difficulties in cross-identifying frequencies detected in 
both seasons. That is why we provide multiple identifications and 
question marks for some of the frequencies in Table [3] 



3.4 Frequencies in the g-mode region 

^From the 2005 data, we recovered all 12 frequencies below 
0.9 mHz found in the combined 2004 data (Table|2j- The frequency 
/g+, an aUas of /g, was detected in the 2005 data and we think 
that it is the true frequency of this mode. In addition, a dozen new 



frequencies in this region were found. We notice that many of the 
amplitudes changed considerably between 2004 and 2005. While 
the 2004 spectrum was dominated by /a and /b with B ampli- 
tudes near 2.7 mma, the 2005 spectrum is dominated by fc with an 
amplitude changing from ~6 mma to ~3 mma (Fig.|4]l. 

The frequencies which we interpret as g-modes occur mainly 
between 0. 1 and 0.4 mHz. However, several periodicities at higher 
frequencies, starting from 0.6 mHz up to /aa at 2.7 mHz, were 
found. Some of them were already detected in the 2004 data. It is 
difficult to prove this, but the comparison with model frequencies 
indicates at least for those with the lowest frequencies, the most 
plausible explanation is also a g-mode pulsation. 
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Table 3. continued. 



# 


Mean freq. 


Amplitude, A [mma 


, dA/dt [mma/d] 




Phase [rad] 






[mHz] 


Uuv 


Bvb 


Vy 


R 


Uuv 


Bvb 


Vy 


R 


/79 


5.468920(22) 




0.26, — 








2.82(15) 






/so 


5.476806(23) 




0.25, — 








4.63(16) 






fsi 


5.477799(25) 




0.24 








2 36n61 






f&2 


5 4S41 97^241 




0.25, 








5.66(15) 






fs3 


5.487352(19) 




0.32, — 








2.58(13) 






/S4 


5.491241(16) 




0.39, 








2.90(11) 






/85 


5.496090(20) 




0.30, — 




0.46, — 




0.77(13) 




1.24(20) 


/se 


5.503194(26) 




0.24, — 








5.63(17) 






/S7 


5.510788(25) 




0.24, — 








4.45(17) 






/ss 


5.516547(19) 




0.30, — 








2.13(13) 






f'i2' 


5.517516(25) 




0.23, — 








0.11(17) 






/89 


5.522630(22) 




Q 28 








s 07n s"! 






/90 


5 531 1 QUI ^'l 




Q 3g 








3 8 1 n 0^ 






/91 


5 534Sn7('99'l 




0.27, 








0.56(15) 






/92 


5.544758(22) 




0.26, 








3.58(15) 






/93 


5.593452(26) 




0.21,— 








3.41(18) 








j.o i'+yzoZ'+ 


A Q7 4_n 095 




1 SQ _i_n DOS 


1 _0 009 


A 09S('9^'l 




4 124(24) 


A 009 ('09 1^ 


A + fo 




1 

1 .oj, 


1 7Q n 01 Q 


1 n HAS 


1 ^0 


5.12(07) 


J .\joy\LZ. ) 






/1 + /3 


5.6323002 


2 99 


9 5S 4-0 01 1 


9 9^ _l_n 091 


9 34 —0 01 3 


fnic\A\ 








/1 + /4 


5.6340636 




0.69, — 


0.89,— 


0.52, — 




0.25(06) 


0.36(11) 


49(17) 


h + h 


5.6478387 




0.23, — 








1.84(17) 






2/3 


5.6496722 




0.33, — 








2.85(12) 






/1 + /6 


5.666189 




0.18,— 








3.65(22) 






3/1 


8.4223924 




0.49, — 




0.60, — 




2.57(08) 




2.63(14) 


2/1 + /2 


8.4379308 




0.21,— 








3.44(19) 






2/1 + h 


8.4397643 




0.22, — 








5.31(18) 








Nobs 


21251 


63497 


17307 


18375 












(TA [mma] 


0.12 


0.04 


0.09 


0.09 












"dA/dt [mma/d] 


0.012 


0.003 


0.010 


0.009 










Detection threshold [mma] 


0.58 


0.20 


0.46 


0.42 










Residual SD [mma] 


12.05 


6.59 


8.58 


8.08 












0.4 0.5 0.6 0.7 
Frequency [mHz] 



1.0 



Figure 5. Fourier spectrum of Bal09 below 1 mHz for B 2004 (top) and 
Bvh 2005 (bottom) data. 



3.5 The region of the dominant frequency (2.8 mHz) 

This is the first of four narrow regions in frequency where p-modes 
occur. The /i (the highest a mplitude term) ha s been previously 
identified as radial mode by ISaran et al.l ( I2OO8I) . The /2, /s and 
/4 form an equidistant triplet, presumably a rotationally split I = 
1 mode (Bar05) and we also recovered /s, /e and in this re- 
gion. Additionally, two new frequencies, /ss and f^j, were found 



E 

£ 1 



< 10- 



2005 
QUINTUPLET 

J 



2004 



2.81 2.82 2.83 2.84 2.85 2.86 

Frequency [mHz] 

Figure 6. Schematic amplitude spectra near 2.8 mHz of 2005 (top) and 
2004 data (bottom). 



in the 2005 data. Together with /s, /g and they form a nearly 
symmetric, although not equally-spaced, quintuplet which is prob- 
ab ly a £ = 2 ro tationally split mode. This was partially confirmed 
bv lBaran et al. 1 12008) using multicolour and spectroscopic obser- 
vations. In addition to the two multiplets and the radial mode, we 
found no other frequencies in the region of 2.8 mHz except the 
'ghost' terms already discussed in Section [373] We also noticed that 
the multiplet spacings changed between 2004 and 2005, becoming 
wider with time. This is shown schematically in Fig. |6] We will 
discuss the possible explanation of this behaviour in Section|5] 

The frequencies of the multiplet components and the splittings 
for 2004 and 2005 are given in Table |4] We can see from this ta- 
ble that while the components of the triplet were within the errors 
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Table 4. Frequencies and splittings for all modes near 2.85 mHz. 







1 ICULICllL/V 


1 icuuciiuy 


Ti 1 TTf^Tf* n r* f 


# 


{e,rn) 


2004 


2005 


in freq. 






[mHz] 


[mHz] 


[mHz] 


h 


(0,0) 








h 


(1,+1) 


2.8232390(06) 


2.8230026(08) 


-0.2364(10) 


h 


(1,0) 




1 S9ilS^(i1 (^OM 

Z..OZ.*tO JU i \UiJ ^ 


4-0 0981 n 

n^w.wz.o 1 i J y 


fA 


(1,-1) 
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equidistant in 2004, this is no longer the case in 2005 and a small 
asymmetry appeared. Moreover, the splittings were larger by about 
15% in 2005. Not only the side components of the triplet changed 
their frequencies, but that of the central one also changed by a sig- 
nificant amount. The changes of frequencies of the triplet compo- 
nents can be even better illustrated in the O— C diagram (Fig.|7) 
where nightly times of maximum light were used. For each fre- 
quency, they were obtained by subtracting all other frequencies 
from the Bvh data and then dividing the residuals into separate 
nights. Then, times of maximum brightness were derived for each 
night. It can be seen that while within the time interval covered 
by observations in each season (40-50 days), the period can be re- 
garded as nearly constant (though different for 2004 and 2005), the 
changes of periods between the two seasons are evident and very 
large for /2 and fi. Unfortunately, it is not known how fast and 
when the periods underwent such a large change. There is an obvi- 
ous ambiguity in counting pulsation cycles between 2004 and 2005 
data. This means that in Fig. [7] the 2005 points could be shifted 
upward or downward by an integer number of pulsation cycles. 

A comparison of the 2004 and 2005 frequencies of the com- 
ponents of the quintuplet indicates that the quintuplet behaved in 
a similar way as the triplet, i.e., increased spacings. The change in 
spacings was different for the |m| = 1 and \m\ = 2 components, but 
as a first approximation we can claim that it was proportional to 
\m\. We may therefore speculate that the same mechanism caused 
the change of splittings in both multiplets. It is interesting to note 
that the main frequency, /i, also changed slightly between 2004 
and 2005. It would be interesting to monitor this behaviour on a 
longer time scale as it could be us ed to determine evolutionary 
changes of the star jReed et al.ll2004h . 

Amplitudes have also changed for these frequencies. The 
largest change was observed for /2 which decreased from 20.4 to 
9.5 mma in BIBvb data, but the amplitudes of /i, /4, /s, /e and /t 
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Figure 7. The O— C diagrams calculated from the nightly times of maxi- 
mum light for the components of the 2.82 mHz triplet. Periods correspond- 
ing to seasonal mean values are given in the panels. 

also decreased from 2004 to 2005. Only for /a did the amplitude 
slightly increase. 



3.6 Other regions with p-modes 

There are three other frequency regions where p-modes occur, 
namely at 3.8 mHz, 4.7 mHz and 5.5 mHz. These regions contain 
frequencies with very small amplitudes. Most of them were de- 
tected only in the Bvb data and very few of these were found in 
the 2004 data. Of the frequencies from 2004 that were recovered 
in the 2005 data, they typically had smaller amplitudes, /g being 
an example. A proper correlation of the 2004 frequencies to those 
in the 2005 data has three complications: (i) The presence of many 
frequencies with similarly small amplitudes; (ii) strong aliasing; 
and (iii) frequencies corresponding to |m| > modes might have 
changed their frequencies considerably between 2004 and 2005, as 
seen in the multiplets. These uncertainties lead to ambiguities in the 
frequency correlation which are indicated in Table [3] using double 
identifications and question marks. 

As the frequency densities are high within each region, some 
of the frequencies are likely members of rotationally split mul- 
tiplets. However, I and m mode assignment would be uncertain 
from our data. Nevertheless, their frequencies may be useful for as- 
teroseismology once the character of long-term frequency-spacing 
changes is established. It is also possible that other, currently un- 
detected components of multiplets will be detected in future ob- 
servations, allowing the full structure of multiplets to be observed. 
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These regions are therefore potentially very important for future 
asteroseismological analyses. 

There is also an isolated frequency (/as) near 3.036 mHz, 
which is between p-mode regions. Its amplitude is very small but 
its detection has been confirmed by G. Fontaine (private communi- 
cation) in an independent set of data, and so we consider it to be 
intrinsic to Bal09. 



3.7 Combination frequencies 

^From the 2005 data, 17 combination terms were detected in five 
distinct regions. They are typically a simple sum or difference of 
two frequencies. Of them, the most important are those that in- 
volve g- and p-modes as it proves that both kinds of modes occur 
in the same star. They were previously detected by Bar05 and are 
detected in the combined 2004 data. There are six g- and p-mode 
combination frequencies observed in the 2005 data and three in the 
2004 data. We do not recover the /s — /i and /i + /g combina- 
tions present in the 2004 data (Table|2j, but we do detect combina- 
tion frequencies in two new regions: sums of g- and p-modes and 
2/i — /a . We also detect two harmonics of /i . Their amplitudes and 
phases indicate that despite a decrease of the B-filter amplitude by 
about 14% between 2004 and 2005, the shape of the light curve of 
the main mode remained practically unchanged. 



3.8 Remarks on tlie frequency spectrum 

Before we try to compare observed frequencies with those de- 
rived from stellar models, a few general comments on the fre- 
quency spectrum need to be made. The substantial data sets of 
2004 and 2005 indicate undoubtedly that amplitude variations are 
very common in this star. Some amplitudes changed by a factor 
of 2-3 over the span of our data. Figure |4] indicates that the time 
scale of amplitude changes might be as short as two weeks. Am- 
plit ude changes were also observed in other sdBV stars. For exam- 
ple, iKilkennv et al.l ( l2006h present the results for four sdBV stars 
which have been observed at least twice. Their analysis showed 
that for three of them the amplitudes of some modes changed sig- 
nificantly. For example, for EC 20338— 1925 the dominant mode 
decreased its amplitude from 26 to 4 mmag. A similar change was 
found for V 338 Ser. Amplitude change s have also been observed 
in other sdBV stars. iReed et al] ( l2007bl) examined follow-up data 
on 23 sdBV stars (including some Bal09 data) and of 54 consis- 
tently detected frequencies, 34 (63%) changed amplitudes by more 
than a factor of two over the duration of their observations (from 
a week up to several years). Fourteen of their frequencies changed 
amplitudes by more than a factor of five; including one that ranged 
from 57 to 8 mma. Therefore, the amplitude variations in Bal09 are 
consistent with those observed in many other sdBV stars. 

We also detect a clear change of frequencies, especially of 
those known to have |m| 7^ 0, i.e., the side components of the ro- 
tationally split multiplets. A small frequency change of the radial 
mode /i was also detected. It would be interesting to trace these 
changes during a longer time scale, not only in Bal09, but also in 
other sdBV stars. They can be used to constrain the evolution of 
sdB stars and/or indicate the mechanism that causes these changes. 
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Figure 8. Physical parameters for Bal 09 derived by different authors. Three 
sdB evolutionai'y sequences (A, B, C) used for theoretical compaiison are 
also shown. He-core burning starts at the point closest to the sequence name 
for each evolutionary track. The lowest gravity phase of a track corresponds 
to the He-core content of ~0. 1 , when the core abruptly engulfs material 
from the radiative zone, causing the loops seen in the tracks at this evolu- 
tionary stage. 



4 COMPARISON WITH THEORETICAL FREQUENCIES 
4.1 Tlie models 

Once the observed frequencies have been determined, it is inter- 
esting to compare them with theoretical frequencies computed for 
models of sdB stars. For this purpose, we computed three sdB evo- 
lutionary sequences (A, B, C). These sequences begin at the onset 
of He-core burning and their parameters are given in Table|5] while 
the evolutionary tracks are shown in the Tch - log g diagram of 
Fig. [8] For comparison, four spectroscopic determinations of T^s 
and log g of Bal09 are also shown with error bars. 

Sequences B and C were obtained by evolving an 1 Mq 
st ar from the main sequence u sing the stellar evolution code 
ofE imenez &MacDonald l l 19961) . A different mass loss Reimer's 
parameter on the red giant branch leads to sdB tracks with differ- 
ent H-envelope mass (A^h), although with very similar total mass 
(~O.47M0), as the core needed to ignite He-core burning is al- 
most the same in all cases. No binary mass transfer was included. 
The position of an evolutionary sequence in the Toff — log g plane 
depends mainly on the total mass (Af) and on A/h. Tracks with the 
same M are located along a diagonal in this plane, ordered by A/h: 
the lower A/h, the larger Tcff and log g are. The diagonal is shifted 
to higher (lower) temperatures if M is increased (reduced). Thus, 
comparing the spectroscopic parameters of Bal09 with the location 
of sequences B and C (Fig. (8)1, we may conclude that either the 
star is at the final stages of He-core burning (or past it) or it is at 
the beginning of the sdB phase, but with a total mass higher than 
0.47 Af0. In the latter case, formation channels involving binary 
interaction should be included to account for the higher mass. To 
allow for the posibility of higher mass for Bal09, a third evolution- 
ary sequence (A) was computed ad hoc with a slightly higher total 
mass (0.55 Mq) than for sequences B and C and was evolved with 
the same code. Its parameters at He-core ignition are also included 
in Tabled 
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Table 5. Parameters of the models at the beginning of the He-core burning 
phase for three evolutionary sequences. Provided are the total mass (Af), H- 
envelope mass (Afjj), effective temperature (T^fi) and logarithm of surface 
gravity (log g). 
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Figure 9. Schematic frequency spectrum of Bal 09 in the region where the 
likely p-modes occur. The dashed line corresponds to /sg . 



4.2 Analysis of p-modes 

In Fig. [9] a schematic spectrum of Bal 09 in the p-mode domain is 
presented. As can be seen, all p-modes except /ag cluster in four 
narrow (from 0.05 to 0.12 mHz wide) frequency groups near 2.8, 
3.8, 4.7 and 5.5 mHz, consisting of 9, 22, 24, and 16 frequencies, 
respectively (Table[3}. 

U sing combined mul ticolour photometry and spectroscopy of 
Bal09. lBaranetal.l j2008h identified /i as a radial mode and con- 
strained the £ value for the 2.825 mHz triplet components and /g. 
For a complete identification of /i , its radial order (n) needs to be 
determined. Assuming that /i is the fundamental radial mode (I 
= 0, n =1) and using the grid of sdB models of ICharpinet et al.l 
( l2002h . Bar05 came to a reasonably good agreement between the 
observed frequency and the location of Bal09 in the log g ~ log T^s 
diagram (see Fig. 13 in Bar05). The first overtone was not regarded 
as an alternative for /i by Bar05, however. Let us consider such 
a case. The period of this fundamental radial mode (Po) would be 
of the order of 356 s/0.74 ~ 480 s, which is not detected in our 
data sets0 Model periods can be matched within the spectroscopic 
log g an d T^g constraints only i f slightly more massive models than 
those o f ICharpinet etal] j2002l) are considered, e.g. some models 
on sequence A in Fig.[8l We can also find an appropriate model if 
we assume that /i is the second radial overtone, though the match 
is worse. In this case Po equals to 593 s. However, all models of a 
given Po have nearly the same value of log g. This can be explained 
as the period of the fundamental radial mode depends primarily on 
average density, i.e., stellar radius. For a set of models that do not 
differ very much in mass this implies that models with the same 
pulsation period would have approximately the same logg. This 
can be seen, for example, in the bottom panel of Fig. 13 of Bar05. 
As the spectroscopic error box for log g and Tcff is still quite large 
(Fig. [8) we cannot constrain n for /i since both the fundamen- 
tal and first overtone are allowable, and even the second overtone 
cannot be excluded. While only higher overtones can be rejected. 



^ 0.74 is the period ratio of the first overtone and fundamental radial modes 
(Pi/Po) found from theoretical models; see also Fig. 1101 



0.9 
0.85 

0.8 
0.75 

0.7 
0.65 

0.6 
0.55 

0.5 



-3.8 mHz O 

e©— i 



Ci«C«!Oo9eSC!D CSS© CSSJO P2/f^ 
— ■ P^/Po 



o«°o9» a»d4Jca88jK»ct®E)«fflog©as c®@) csso Pj/F^ 




-5.5 mHz ^ 



100 



200 



300 



400 



500 600 

Po[s] 



700 



800 900 



Figure 10. Period ratios for the first four radial orders of radial modes plot- 
ted as a function of the period of the fundamental radial mode, Po. Two 
sets of models were used: (i) those described in Section l4~nand listed in Ta- 
ble|5](large symbols), (ii) the sequences 5 to 7 from Ch arpinet et alj ^20021) 
(small symbols). Three horizontal strips denote the observed ranges of fre- 
quency ratios. All modes from p-mode groups near 3.8, 4.7 and 5.5 mHz 
are considered. 



we can constrain \ogg for our three possibilities to 5.51 (/i fun- 
damental), 5.34 (1st overtone), or 5.22 (2nd overtone). Improved 
constraints on log g will distinguish between these three possibili- 
ties. However, we prefer to assume that /i is the fundamental radial 
pulsation as this is the most likely possibility. If /i were the first 
overtone then the fundamental radial pulsation (and also some non- 
radial p modes with similar frequencies) should occur near 2. 1 mHz 
and these pulsations are not observed. 

Another argument in favour of /i as the fundamental ra- 
dial mode is the overtone ratio. If we assume that at least one 
mode in the 3.8 mHz group is also radial then the overtone ratio 
is 2.807 mHz/(3.774^ 3.841) mHz « 0.731-^0.744 which can be 
compared with theoretical ratios, as is done in Fig.[TO] It can be seen 
that for most models the Pi /Po ratio is near to 0.74 and only for 
model C is also P3 / Pi allowed in some of its evolutionary phases. 
The latter case implies fi to be the first overtone. As such, we can 
conclude that if any of the 3.8 mHz-region frequencies are radial, 
then /i is most likely the fundamental radial pulsation. Addition- 
ally, Fig.[TO]indicates that by choosing fi to be the fundamental ra- 
dial mode, the 4.7 mHz and 5.5 mHz groups match overtone ratios 
for P2 / Po and P3/P0. It seems unlikely that such a good match is 
serendipitous. Reversing this reasoning we may therefore conclude 
that each of the groups of p-modes in Bal09 may contain a radial 
mode (if so, we know which overtone it is). Alternatively, the same 
could be true for consecutive overtones of non-radial modes as they 
have period ratios similar to the radial ones, imless when avoided 
crossing occurs. 

Considering splittings, we have searched for additional multi- 
plets using the same frequency spacing as observed for multiplets 
at 2.8 mHz. Because of the high frequency density in these regions, 
there are likely rotationally split multiplets. However, the simple 
structure of frequencies in the 2.8 mHz group contrasts with more 
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Figure 11. A schematic frequency pattern in four groups of p modes: near 
2.8, 3.8, 4.7, and 5.5 mHz. 



complicated frequency patterns in the three remaining groups of 
p-mode frequencies (Fig. lilt. Another hkely multiplet is the sym- 
metric triplet of fr2, /ts, and fj^ (TableO with splittings of 1.440 
± 0.021 and 1.446 ± 0.023 /iHz. Bearing in mind the richness 
of the spectrum, observed amplitude variations, and uncertainties 
caused by aliasing, we prefer not to speculate further on additional 
multiplets until better data more fully reveal the frequency spec- 
trum. The high frequency density also indicates that all groups of 
p-modes must have £ > 2 modes, except those in the 2.8 mHz re- 
gion. 



4.3 Analysis of (?-modes 

Bal09 shows a very rich spectrum in the low-frequency do- 
main, which theoretically corresponds to (j-modes. In Fig. [T2l the 
schematic amplitude spectrum of Bal 09 in this frequency range is 
displayed. For comparison, theoretical frequencies of a representa- 
tive sdB model in sequence C are also shown. We did not intend to 
find a model which matches the observed frequencies, but just to 
compare the density of theoretical and observed frequencies in this 
region. As can be seen in Fig. [12] at least £ = I and 2 modes are 
needed to account for the observed frequencies. However, as the 
density of theoretical frequencies increases with £, the observed 
spectrum would be better reproduced if higher degree modes were 
considered. 

We also notice a lack of frequencies between 0.4 and 0.6 mHz. 
Either they have amplitudes below our detection threshold or they 
are not excited. We investigate the possibility that they are not 
excited and show in Fig. [T3] the derivative of the work integral 
(dW/dr) for g-modes with £ = 2 and radial orders \n\ = 1,7 and 
30. These modes were chosen as representative of low, intermedi- 
ate and high radial orders. We used the Gra Co non-adiabatic code 
jMova et al.l |2004| : iMova & Gairidol l2008h to compute dW/dr, 
which is positive in regions contributing to driving oscillations. 
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Figure 12. Schematic frequency spectrum of Bal 09 below 1 mHz com- 
pared with theoretical ones for g-modes with £ = 1,2 and 3. The latter are 
shown for a representative sdB model from sequence C. Numbers on the 
right indicate the value of n for the first g-mode below 1 mHz. 
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Figure 13. Derivative of the work integral for £ = 2, n = —I, —7, and 
—30 g-modes as a function of a fractional stellar radius for a representative 
model from sequence C. Regions contributing to driving (damping) produce 
positive (negative) dW/dr values. The Rosseland mean opacity is included 
as a solid line and plotted in arbitrary units. The gray area delimits the 
driving region associated with an iron-group opacity bump. 



while negative in damping regions. For hot B subdwarfs the driving 
zone is related to the iron-peak elements opacity bump. In Fig. 1131 
the Rosseland mean opacity for theoretical model C used is plot- 
ted as a solid line. Fig. [T3] shows also that, while the region of 
the bump (shown in gray) contributes to the excitation of low- and 
high-radial order (;-modes (dW/dr > there), it has a negligi- 
ble effect in driving intermediate-radial order modes. This result 
could support the hypothesis that g-modes with intermediate radial 
orders are not excite d in Bal 09. Non-adiabatic computations by 
Ijefferv & Said ^200^ also result in exciting low- and high-order g- 
modes but damping intermediate ones for some models. This seems 
a reasonable explanation for the lack of observed frequencies be- 
tween 0.4 and 0.6 mHz. Theoretical g-modes with low radial orders 
can have frequencies up to ~2.3mHz (not shown in Fig.ll2t. We 
detect three frequencies {fx, /y and /z) between 1.8 and 2.3 mHz 
which could then be attributed to low radial order g-modes. 

Most g-modes observed have periods in the range between 
2500 and 7200 s, although a few extend this range to 1200-10300 s. 
Theoretical (?-modes with £ = 1-3 and periods in the range 2500- 
7200 s (Fig. 112b have radial orders in the range |n| ~ 1 1-80, which 
extends to |n| ~ = 5-115, if we include all of the observed peri- 
ods. With high-radial orders, we can ex pect asymptoti c behaviour 
which results in eq ually-spaced periods (Ffassoullll98Q) . However, 
iKawaleJ ( Il999h and lCharpinet et al.l ( 12000 j) have shown that sdBs 
might display trapped g-modes, mainly caused by the chemical 
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Figure 14. Schematic g-mode spectrum in periods of Bal 09. Period dis- 
tances between adjacent modes are indicated with numbers expressed in 
seconds. 

transition between the He and H radiative layers in the envelope. 
Mode trapping causes some periods to move away from asymp- 
totic behaviour. In Fig. [14] we show a schematic g-mode spectrum 
in periods and indicate the differences (in seconds) between adja- 
cent periods. However, the observed period spectrum does not ap- 
pear equally spaced, which suggests that some n orders were not 
detected or excited. Fig. [12] also shows that the density of theoret- 
ical periods is much higher than the observed density, even if only 
periods with ^ ^ 3 are considered. 

We searched for equally-spaced periods using different tech- 
niques. While we could produce some reasonable alignments in 
echelle diagrams (e.g. for period spacing AP ~ 286 s), the results 
are not convincing. If mode trapping is effective, we should not ex- 
pect good alignment from an echelle diagram. A way around this 
problem would be to fine-tune models using the p-modes and then 
search the g-modes for periodicities, or even match their periods 
directly to the model ones. If some p-modes could be identified, it 
could greatly increase our understanding of sdB interiors. 



5 DISCUSSION AND CONCLUSIONS 

The 2005 campaign on Bal09 was carried out because the object 
is one of the most interesting pulsating sdB stars. It has high- 
amplitude oscillations, it is one of the brightest sdB pulsators {B 
= 11.8 mag), and it has the richest pulsation spectrum of both p- 
and (?-modes. The 2005 campaign led to the discovery of nearly a 
hundred independent frequencies (TableJSj. Moreover, we showed 
that, like in many other sdBV stars ( iReed et al. 2007b), amplitudes 
of many modes in Bal09 vary and the time scale of these changes is 
as short as several weeks. The 2005 data also allowed the detection 
of all components of the quintuplet. This was possible because of 
the low detection threshold of ~0.2mma in the Bvb data, which 
also resulted in the detection of a large number of frequencies, es- 
pecially near 3.8 and 4.7 mHz. Even if some of them turn out to 
be non-axisymmetric multiplets, the frequency density is such that 
at least some of them must have I > 2. The most surprising result 
of this campaign is the discovery of a considerable change of split- 
tings in two multiplets (a triplet and quintuplet) which are believed 
to be rotationally split 1= \ and 2 modes. As far as we are aware, 
this is the first clear detection of such a change of splitting in any 
pulsating star. 

The presence of multiplets is usually attributed to rotation of 



the star, causing the m degeneracy to be broken. Two clear mul- 
tiplets are detected in Bal 09: a £ = I triplet at 2.82 mHz and a I 
= 2 quintuplet at 2.85 mHz. The structure of the quintuplet is es- 
pecially interesting. Its components, all of which are detected in 
2005, form a symmetric but not equidistant pattern. Specifically, the 
mean separation between |m| = 2 and |m| = 1 components (2.100 
± 0.014 fiHz) is considerably larger than the mean separation be- 
tween |mj = 1 components and the central peak of the quintuplet 
(1.577 ± 0.007 fiHz). The same was true in 2004. Such a pattern 
cannot be produced by a solid body rotation, nor even by a spher- 
ically symmetric rotation, f2 = f2(r). Indeed, in both these cases 
splitting within a multiplet is, in the limit of slow rotation, propor- 
tional to m (e.g. iKawaler & Hostletll2005h . Consequently, compo- 
nents of the multiplet have to be equidistant. This is not what we 
observe. Thus, the observed structure of ^ = 2 quintuplet implies 
that rotation rate of Bal 09 must depend on the stellar latitude, 6. 

Additional information can be inferred from comparison of ro- 
tational splitting for |m| = 2 and for |mj = 1. Because of its spacial 
structure, modes of |m| = 2 are most sensitive to equatorial rota- 
tion of the star, while modes of |m| = 1 sense mostly rotation in the 
intermediate latitudes. The observed splitting is larger for |m| = 2, 
which implies that rotation of Bal 09 is fastest on the equator. This 
has to be true at least in a large part of the star's envelope. Quali- 
tativ ely, this dependence on 9 is similar to that inferred for the Sun 
(e.g. lThompson et al.ll2003h . 

The most puzzling result of 2005 campaign is the detection of 
changing rotational splittings. Both in the triplet and in the quin- 
tuplet the separation of components is much wider in 2005 than in 
2004. The differences range from 12% for ^ = 2 |m| = 1 to 16% 
for 1 = 2 \m\ =2. Similar behaviour of both multiplets puts strong 
constraint on any theoretical explanation of the observed effect. For 
example, it excludes internal 1:1 resonances within the multiplets 
( iBuchleret ai]|l995h . Such a coupling could modulate the separa- 
tion of multiplet components, but it would act in each multiplet in- 
dependently, giving no explanation why they all change in the same 
way. Simultaneous widening of both multiplets must be caused by 
a common mechanism, which affects the structure of the star. 

Taken at the face value, widening of multiplets in Bal 09 
would imply a considerable change of its internal rotation, occur- 
ring in just one year. At first glance this seems to be a very far 
fetched hypothesis. However, variations of internal stellar rotation 
on time scale of a few years are not un heard of. Such variations are 
actually observed in the Sun ( IShibahas hi 2004; Howe 2008). These 
so-called torsional oscillations cause periodic changes of the solar 
rotation pattern with the 1 1-year solar cycle. 

The cause of torsional oscillations in the Sun i s not clear, 
although several hypotheses have been put forward jShibahashil 
[2004). We must note that the internal structure of Bal 09 and the 
Sun is very different. In contrast to the Sun, the envelope of an 
sdB star is almost entirely radiative, with the exception of an ex- 
treme ly thin convective laye r associated with the Hell ionisation 
zone i lChan3inet et al.l2000b l). However, the presence of convection 
is not necessary to drive torsional oscillations in the star, as they can 
be ex cited by interaction between dif ferential rotation and magnetic 
field iGoode & Dziembowski|[T99lh . In this mechanism, differen- 
tial rotation and the toroidal component of the field exchange en- 
ergy, causing periodic redistribution of angular momentum inside 
the star. This in tur n, leads to periodic change s of rotational split- 
tings. According to lKawaler & Hostler! ( |2005|) . strong differential 
rotation is what we expect in envelope of every sdB star. The mag- 
netic field in Bal 09 has not yet been measured, bu t observations 
of several other sdB stars reveal fields of '~1.5kG ( lO'Toole et al.l 
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I2005I) . The amplitude of torsional oscillations discovered in the Sun 
is rather small, of the order of ±1 - 1.5% of the local rotation ve- 
locity. This is amount too small to account for changes in multiplet 
splittings as large as those observed in Bal 09. Torsional oscillations 
in Bal 09 have to be at least an order of magnitude stronger than in 
the Sun. Whether torsional oscillations in Bal 09 can be driven to 
such amplitudes remains an open question until detailed model cal- 
culations become available. 

These results indicate that Bal09 may be a key object in our 
understanding of pulsations in sdBV stars and its variability is def- 
initely worth further study. First, the changes of splittings should 
be monitored on a longer time scale as this may help to under- 
stand their origin. The same is true for amplitude changes. Next, 
as the amplitudes change so distinctly, it may happen that modes 
undetectable in one season will become visible in the other, thus 
allowing to a complete spectrum of excited modes. Finally, the ob- 
served frequencies should be compared to a larger grid of mod- 
els using id entifications of the strongest frequencies as done in 
iBaran et al.l ( I2OO8.I . Such constrained modes should result in con- 
straining global stellar parameters for Bal09 and reveal its internal 
structure. 
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